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Abstract 

We study the problem of homogenization for inertial particles moving 
in a periodic velocity field, and subject to molecular diffusion. We show 
that, under appropriate assumptions on the velocity field, the large scale, 
long time behavior of the inertial particles is governed by an effective dif- 
fusion equation for the position variable alone. To achieve this we use a 
formal multiple scale expansion in the scale parameter. This expansion 
relies on the hypo-ellipticity of the underlying diffusion. An expression for 
the diffusivity tensor is found and various of its properties studied. In par- 
ticular, an expansion in terms of the non-dimensional particle relaxation 
time r (the Stokes number) is shown to co-incide with the known result 
for passive (non-inertial) tracers in the singular limit r — > 0. This requires 
the solution of a singular perturbation problem, achieved by means of a 
formal multiple scales expansion in r. Incompressible and potential fields 
are studied, as well as fields which are neither, and theoretical findings 
are supported by numerical simulations. 



1 Introduction 

Understanding the transport properties of particles moving in fluid flows and 
subject to mole cular diffusion is a problem of great theoretical and practical 
interest [111 123) . For the purposes of mathematical analysis, the fluid velocity 
is assumed to have some statistical or geometrical structure which mimics fea- 
tures of real fluid flows, and yet is such that the resulting equation describing 
the motion of the particle is amenable to both analysis and efficient numerical 
investigations. For example the velocity may be assumed steady and periodic 
in space, or may be assumed to be a Ga ussi an random field in space-time. This 
is the problem of turbulent diffusion [7| ■ 

It is often the case that the velocity field of interest is active at various length 
and time scales. Consequently, the equations which govern the particle motion 
are very hard to analyze directly. In such cases an effective equation which 
governs the behavior of the particles at long times and large scales compared to 
those of the fluid velocity is sought. The derivation of such an effective equation 
is based on multiscale/homogenization techniques 0]. 
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This problem has been studied extensively over the last thirty years for 
passive tracers, i.e. massless particles. It has been shown that, for periodic 
or random velocity fields with short range correlations, the particles perform 
an effective Brownian motion. The covariance matrix of this Brownian motion 
- the effective diffusivity - is computed through the solution of an auxiliary 
equation, the cell problem. In the case of steady periodic velocity fields the cell 
problem is a linear elliptic PDE with periodic coefficients. Various properties 
of the effective diffusivity have been investigated. In particular, it has been 
shown that the diffusivi ty is alway s enhanced (over bare molecular diffusion) 
for incompressible flows [121 1231 and always depleted for potential flows |46j . 
Extensions of the above results to the case where the molecular diffusivity is 
modelled as colored noise have also been analyzed 

However, there are various applications where the particles cannot be mod- 
elled as tracers and inertial effects have to be taken into account. The need for 
adequate modelling and analysis of the motion of inertial particles in various 
applications in science and engineering has been recogniz ed over the last few 
years with a number of publications in this direction |SJ As well as 

bringing important physical features into the problem, the presence of inertia 
also renders the mathematical study of the resulting equations more delicate. 
The purpose of the present paper is to study the long time, large scale behavior 
of inertial particles moving in periodic velocity fields, under the influence of 
molecular diffusion. Both analytical and numerical techniques are used. From a 
physical standpoint the main interest in this work is perhaps the enormous en- 
hancement in diffusion which can be introduced through relatively small inertial 
effects; see section |SJ From the mathematical standpoint, the main interest is 
perhaps (i) the hypo-elliptic structure of the differential operators arising in the 
multi-scale expansions (section [3} and (ii) the elucidation of the passive tracer 
limit, from within the inertial particles framework, again by use of multiple 
scales expansions (section EJ- 

We consider the diffusive behaviour of trajectories satisfying the equation 

tx = v(x) — x + a$, (1) 

where (3(t) is a standard Brownian motion in R d , d > 1. In Figure Q we plot a 
sample trajectory of Q with dimension d = 2, r = 1 and a — 1; the velocity 
v(x) is given by the Taylor-Green field: 

v(x) = V^VtgOe), ^tg(x) = sin(a;i) sin^)- (2) 

This figure strongly suggests that the long time behavior of the particles is 
diffusive. Our goal is to make this type of observation precise. We provide a 
mathematical formula for the effective diffusivity, and study its properties by 
a combination of analysis and Monte Carlo simulations. The results raise a 
number of interesting questions of physical interest, regarding behaviour of the 
effective diffusivity, and also provide a framework for mathematical analysis to 
further probe questions of physical interest. 

Problems of this type have already been studied in a variety of situations: 

1. When a = 0, V • U{x) = and 

v(x) = U{x) + 5tU{x) ■ VU(x) (3) 

we obtain the following model fo r th e motion of a particle in an incom- 
pressible steady velocity field [271 : 

x = 5U(x)-VU{x) + -{U{x)-x). (4) 
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Here x(t) denotes the particle position and U(x) is the fluid velocity. The 
parameter S — — is the ratio between the fluid density and the particle 

density and r denotes the Stokes number. Equation J3J becomes Stokes' 
law when (5 = 0, i.e. when the particles are much denser than the surround- 
ing fluid. The model Q was analyzed numerically for two dimensional 
steady incompressible cellular flows by Crisanti et al. [o| . Their numerical 
experiments showed that for particles slightly denser than the surround- 
ing fluid (1 — S w and S < 1), and for t = 1, the particles perform 
an effective Brownian motion with a well defined diffusion coefficient The 
term inertial diffusion was introduced to describe this phenomenon. In 
Figure 2 we present a sample trajectory of @ with U(x) given by the 
Taylor Green velocity field for r = 2.0 and 5 = 0.7. The figure again 
suggests long time behavior of the particle trajectories which is diffusive. 
This phenomenon has been observ ed and analyzed-mostly numerically-by 
various authors [HI 1261 13*71 1471 148) . The range of 6 in which Brownian mo- 
tion behaviour is observed appears to be connected to the linear stability 
properties of the equilibria of (J3J; see the discussion in section H3 

2. When a = 0, V • U(x) = and v(x) = U(x) (i.e. S = in Q) a similar 
investigation was carried out by Wang et al. in 07] in the case d = 3. 
They studied equation numerically for the three-dimensional ABC 
flow. It is well known that the streamlines of this flow are chaotic 0. It 
was exhibited numerically that, for sufficiently small r, the particles per- 
form Brownian motion. The resulting diffusion coefficient was evaluated 
numerically and its dependence on r was investigated. 

3. When er > and v(x) — —\7V(x) it is possible to study the effective 
behavior of particles moving in a periodic or random potential flow, under 
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Figure 2: The two components of the solution of l@J, for U(x) given by (J2J), 
with S = 0.7 and r = 2.0. 



the influenc e of molecular diffusion. This has been addressed by various 
researchers |15M22ll3"Tll35| . In this case the particle motion is governed by 
the Langevin equation: 



where (3(t) is the standard Brownian motion in R d , a > and V(x) is a 
potential which can be eith er periodic or random. It was shown by Papan- 
icolaou and Varadhan that, for sufficiently smooth random potentials 
V(x), the particles perform an effective Brownian motion, at long length 
and time scales, with a nonnegative-definite effec tive diffusivity. A similar 
result was prove d for pe riodic potentials in |35| and various refinements 
were analyzed in 15 , 22 . The analysis of equation JSJ is greatly simplified 
by the fact that the explicit form of the invariant distribution associated 
to the particle position and velocity is known. 

4. The problem of averaging for equati ons of the form with periodic 
coefficients was studied by Freidlin in |T3]- He considered the equation 



where v(x), u(x) are smooth, 1-periodic functions and the matrix a(x) — 
a(x)a(x)* is non-degenerate. It was shown in J3j that as r, e tend to 
0, the particle position x converges to the solution of a first order SDE 
with constant coefficients, leading to linear in time mean flow, with super- 
imposed Brownian fluctuations. The averaged coefficients of the limiting 
SDE depend on how fast e tends to relative to r and two different cases 
have to be distinguished. 



x = — 



VV{x) -x + a/3, 



(5) 
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The purpose of the present paper is to study the long time, large scale behavior 
of inertial particles satisfying when the velocity field is periodic with period 
1: v(x + ij) = v(x), j = 1, . . . , d where {ej} d =1 are the unit vectors in E d . We 
show, by applying multiscale techniques to the backward Kolmogorov equation 
associated to QJ, that, provided the drift term is centered with respect to the 
invariant distribution of the stochastic process corresponding to Q, the long 
time, large scale behavior of the inertial particles is governed by an effective 
Brownian motion. The effective diffusivity tensor K, is computed through the 
solution of the cell problem which is an equation posed on M. d x T d , where T d 
is the d-dimensional unit torus. It is shown that this tensor is nonncgative and 
that, consequently, the effective dynamics is well-posed. We emphasize that our 
analysis does not assume any specific structure on the drift term and that, in 
particular, it is valid for both divergence free as well as potential flows, provided 
that a centering condition on v(x) is satisfied. However, in the case of potential 
flows, the effective diffusivity is shown to be depleted, for all Stokes' numbers, 
when compared with bare molecular diffusion; this generalizes a known result 
for passive tracers. 

Our theoretical findings are supported by numerical simulations which probe 
the dependence of the effective diffusivity on parameters such as r and, when 
the velocity field is given by , 8. The equations of motion are solved for 
the Taylor-Green flow (J2J) and the long time behavior of the particle trajecto- 
ries is shown numerically to be that of a Brownian motion. The numerically 
computed effective diffusivity is compared with the enhancement for passive 
tracers, i.e. when r = 0, and it is shown that the presence of inertia enhances 
the diffusivity beyond the enhancement for the passive tracers. The problem is 
also analyzed through a small r expansion of the effective diffusivity. We show 
that, to leading order in r, the effective diffusivity is equal to that arising from 
the homogenization of passive tracers. We also compute the 0(r) correction to 
the effective diffusivity for one dimensional gradient flows of the type JSJ and 
show that it is always negative. Thus, in the one dimensional case and for r 
sufficiently small but positive, the diffusivity is depleted over bare molecular 
diffusion even further than the depletion which occurs when r = 0. 

Here we analyze the problem of homogenization of inertial particles and 
derive the formula for the effective diffusivity using formal multi-scale calcula- 
tions. We emphasize that the basic homogenization result derived in this paper 
can be proved rigorously using techniques from stochastic analysis, in particular 
the martingale central limit theorem (l8| : we refer to for the rigorous proof 
of the homogenization theorem. 

The paper is organized as follows. In section [3 we review homogenization 
for passive tracers. The derivation of the homogenized equation for the general 
second order stochastic differential equation with periodic velocity field v(x) 
is discussed in section[3] Various properties of the effective diffusivity are derived 
and analyzed in section 01 The homogenization result for the equation Q , 
subject to molecular diffusion, is also presented there. The small r expansion 
for K- is studied in sectional Numerical experiments for the Taylor-Green flow 
are presented in section Section [7] contains our conclusions. Finally, some 
technical results which are needed for the rigorous justification of the multi-scale 
method employed in this paper are proved in the appendix. 



2 Homogenization for Passive Tracers 

In this section we review the homogenization result for passive tracers - i.e. 
massless particles - advected by a velocity field v(x) and subject to molecular 
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diffusion: 

x = v(x) + a$. (6) 

We do this in order to emphasize the structural similarities between this case 
and the inertial case QJ which is the focus of the paper. 

The velocity field v(x) is assumed smooth and periodic with p eriod 1. We 
will take a > 0. This problem has been extensively analyzed, cf. [41 1281 l2~9l 14fi| . 
Here we merely outline the main steps in the derivation of the homogenized 
equation. Notice that we do not assume that the flow is either incompressible 
or potential. We also remark that we have chosen to base our analysis on 
the Backward Kolmogorov equation as opposed to Forward Kolmogorov-i.e. 
the Fokker-Planck- equation, which is more customary. Of course, the two 
approaches are equivalent, provided that the field v(x) is sufficiently smooth. 
We feel, however, that the analysis based on the backward Kolmogorov equation 
is cleaner than that based on the Fokker-Planck equation. Furthermore, the 
backward Kolmogorov equation i s the st arting point for the rigorous justification 
of the results reported here, see [THIITHj . 

The process x(t) whose evolution is governed by JfjJ) is ergodic on the torus 
T d , ch. 3] for details. This in particular implies the existence of a unique, 
smooth, invariant density ttq(x) which is the unique solution to the stationary 
Fokker-Planck equation. We will assume that the drift term averages to with 
respect to this invariant density ttq: 

v(x)ir (x) dx = 0. (7) 

The above centering condition ensures that the long time, large scale behavior 
of the particle is diffusive and, in particular, excludes the possibility of the 
existence of an effective drift. It is automatically satisfied for potential flows 
(for which the invariant density ttq(x) is the Boltzmann distribution) and it 
reduces to the condition that the average of the drift over the unit cell vanishes 
for incompressible flows, since in this case ttq(x) = 1 1 . 

2.1 The Rescaled Process 

We use the scaling property of Brownian motion f3(ct) — y/c/3(t) in law |171 
Lemma 2.9.4] to deduce that, under the re-scaling t — ► t/e 2 , x — > x/e we obtain 

1 x 

X = -v(-) + (7/3. 

Setting z = x/e gives the system of SDEs 

x = -v(z) + cr/3, (8a) 

z = ^v(z) + ^/3, (8b) 

1 It would seem more natural from a physical point of view to impose the condition 
Jjd v(x) dx = 0, rather than condition J7J, since the invariant density itself depends on the 
drift v(x). But, as we will see, condition Q is natural from a mathematical point of view. 
If, instead, the condition J Td v(x) dx = is imposed, then an effective drift can appear and 
then a Galilean transformation with respect to this effective drift can bring the problem to 
the form that we consider. 
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with the understanding that zeT and x 6 Mr. This clearly exhibits the fact 
that the problem possesses three time scales of 0(1), 0(e) and 0(e 2 )- We now 
average out the fastest scale, given by z. 



2.2 Mutiscale Expansion 

The backward Kolmogorov equation for J5J) is 



' 1 ;C Q + -A + £2 ) u, (9) 



c 



dt 
with 

2 

C Q = w(z)-V z + yA z , 
£1 = v(z) ■ V x + a 2 V z .V x , 
C 2 = yA x . 

Here Co is a uniformly elliptic operator with periodic boundary conditions on 
T d . We look for a solution of equation which has the form: 

u = u + eu\ + e 2 u 2 + . . . 

with m = Ui(x, z,t),i — 1,2,.... Substituting the above ansatz into we 
obtain the following sequence of equations: 

£ou = 0, 

Cqux = -£i« , (10) 

£ou 2 = -C\ux -£ 2 u + -5—. 

at 

The ergodicity of the process generated by Cq implies that the null space of the 
operator Cq consists of functions which are constants in z. We let ttq be the 
invariant density which satisfies the equation 

/> (z)=0. (11) 

Thus the equation Cof = g has a solution if and only if g averages to over T d 
with respect to 71-0(2:): 



(3)770 : = / g(z)n a (z)dz = 0. (12) 

We refer to ch. 3] for proofs of these results. 

In view of the discussion in the previous paragraph, from the first equation 
in (|10f) we deduce that uo — uo(x,t). Furthermore, the centering condition on 
the drift term v(x) implies that 

C 1 u (x,t)TT (z) dz = I / v(z)n (z) dz I • W x u = 0, 

T d \JT d / 
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and consequently the second equation in llUfl is well-posed. We can solve this 
equation using separation of variables: 

Ui(x,z,t) = x(z) ■ V x Uo(x,t) 

where the corrector field x( z ) satisfies the cell problem : 

Cox(z) = -v{z). (13) 

Notice that, in view of the centering condition on the drift term, the cell problem 
is well-posed and it admits a unique solution. 

Now we proceed with the final equation in (|10fl . We apply the solvability 
condition to obtain: 

= J ^CiUx + C 2 Uo - ~7^j tto(z) dz 

Oil £7^ ( f \ 

= —j£ + yA.no + (j [v(z) ® x(z) + a 2 V zX {z)} tt (z) dz) : D 2 x u 
= -—+K:D x u . 

We have used the notation D 2 u for the Hessian of Uq: D 2 uq = ( dx -fa ) 

and A : B for the product of the matrices A and B. Moreover, ® stands for the 
tensor product between two vectors. 

From the above equation we deduce that: 

8t K, 'i).r,i).r/ 1 ' 

where the summation convention has been used. The effective diffusivity K is 

a 2 f f dxi 

£ij = ~TT S v +/ Vj(z)xi(z)ir (z)dz + a / — -(z)-kq(z) dz. (15) 

- JT d JT d 

Equation l|14|) is the backward Kolmogorov equation associated to a Brownian 
motion X(t) with covariance matrix /C. We remark that these formal calcula- 
tions can be justified rigorously using either energy esti ma tes 0], the method 
of two-scale convergence pQ or probabilistic methods E2 • 

2.3 Well-posedness of the Limiting Equation 

The limiting backward Kolmogorov equation (|14|l is well-posed, i.e. the effective 
diffusivity is a nonnegative matrix. To see this, we first observe that 

W*o) = -7T C f + <J 2 (A z f)7r Q +a 2 V z f- V z tt 

= -vroA,/ + <t 2 V 2 -{(V,/)7t }, (16) 
for every smooth periodic function f(z). Consequently: 

/ /(£>/) 7r <fe = -2- f |V 2 /| 2 vr dz, (17) 

Jfd I Jfd 
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in view of (|16[) and an integration by parts. Now let a be an arbitrary vector 
in R d and let \o — a -X where \ is the solution of the cell problem l|13|) . The 
scalar quantity xo satisfies the equation 

£0X0 = -a -v. (18) 

We consider the effective diffusivity along the direction a. We use (|17|) to obtain: 

a-Ka = ykl 2 + / T d(a ■ v)xa^odz + a 2 f Td (V z Xo ■ a)ir dz 

= j- \a\ 2 - / T d {£oXo)xoKo dz + cr 2 J Td (V z % • a)ir dz ^ 
= ^-\a\ 2 + ^ J Td \V z Xo\ 2 Tr Q dz + (T 2 J Td (V z xo ■ a)n dz 
= T" It" I V zXo + a\ 2 ir a dz. 

Using the notation i|12[l . the penultimate line in 11911 may be written as 
2 2 

a . fc a = ^_| a |2 + l-(|V,xo| 2 ),* + a 2 (V zX o • o)» . (20) 

From l|19f) we deduce that the effective diffusivity is indeed nonnegative and the 
well posedness of the effective equation ()14|l is demonstrated. 



2.4 Incompressible Flows 

Whether the diffusivity is enhanced or depleted depends on the specific prop- 
erties of the periodic drift term. For the case where the flow is steady and 
either divergence free or potential more detailed information can be obtained. 
In particular, for incompressible flows we have that tto(z) = 1. Con sequently 
the last integral on the right hand side of the penultimate line in 1|19|) vanishes 
on account of the periodicity of xo(z). Thus, the effective diffusivity along the 
vector a becomes: 

2 2 r 

aXa = ^-\a\ 2 + G — / \V zX o\ 2 M z ) dz. (21) 

z l Jj d 

This shows that tran spor t is always enhanced over bare molecular diffusion, for 
incompressible flows |29| . 



2.5 Potential Flows 

When v(z) — —W(z), the invariant density tto(z) is the Boltzmann distribu- 
tion: 

12 f 2 

7To(z) = ^exp( -F(z)), Z= exp( -V(z))dz. 

Z <7 Z Jjd (T z 

An integration by parts, together with the periodicity of Xo( z ), 7r o(^) and equa- 



ntcg r 
tion {T7J, gives: 



\ ^zXo(z).air (z)dz = -a / \V z Xo(z)\ ir (z) dz, 



and consequently, from the penultimate line in (|19J) . we find that 

(22) 



aXa = ^\a\ 2 - ^-/ Td \\7 zX o\ 2 TTo{z)dz, 
^\a\ 2 -^(\V zX o\ 2 ), a 
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with ttq(z) being the Boltzmann distribution. This shows that transport is 
alw ays depleted, compared with bare molecular diffusion, for potential flows 

B3- 



3 Homogenization for Inert ial Particles 

In this section we will derive the homogenized equation which describes the 
motion of inertial particles at long times and large scales using multi-scale tech- 
niques. The equation of motion for the inertial particles is JTJ. 



3.1 The Rescaled Process 

We start by performing a diffusive rescaling to the equations of motion JTJ: 
t — » t/e 2 , x — > -. Using the fact that f3(ct) = y/c(3{t) in law we obtain: 

re x = -v(-) — x + op. 

Introducing y = ^/rex and z — x/e we write this equation as a first order 
system: 

± = ir^ 

y = -jh< z )-T?y + ^i (23) 

i = 

with the understanding that z £ T d and x, y £ M. d . This clearly exhibits the fact 
that the problem possesses two time scales of 0(e) and 0(e 2 ). We now average 
out the fastest scales, on which (y, z), evolve, and show that the fast and large 
fluctuations in x induces diffusion on time-scales of 0(1). 



3.2 Multiscale Expansion 

The backward Kolmogorov equation associated to equations i|23|) is 
du e 1 _ c 1 / 1 _ 1 , s _ 1 



dt yfFe' 



where: 



T" " \/T ' ' T 



y ■ V x u e + - —=y ■ V z + - } =v{z) ■ V y + -C L 



1 Co + -£i ) u\ (24) 



e 2 e 



C = ^ r (yW z + v(z)-W v ) + -C ou , 

VT T 
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Note that C ou is the generator of a standard d-dimensional Ornstein-Uhlenbeck 
process [141 ch. 3]. This process is ergodic with Gaussian invariant density 
satisfying 

(C ou )*p ou = 0. (25) 

In order to carry out the analysis which follows we will make use of the ergodic 
properties of the solution to (Q with x — > z. Using the tools developed in 
one can prove that the process z(t), y(t) with y(t) — y / ri(i) is ergodic on 
T d x K d . 2 The analysis implies that there exists a unique invariant density 
p(y, z) with support of positive measure on T d x R d . The hypo-ellipticity of Cq 
established in 25 shows that the density is smooth and is hence the unique 
solution to the stationary Fokker-Planck equation associated to the process (pQl 
|2H ch. 11]: 

C* p(y, z):=-^=.{y V zP + v{z) ■ V yP ) + i (V y • (yp) + yA 9 j p = 0. (26) 

The Fokker-Planck operator Cq is the adjoint of the generator of the process 
£ Q . The null space of the generator C consists of constants in z, y. Moreover, 
the equation Cof — g, has a unique (up to constants) solution if and only if 



(g)p : = / g(y,z)p{y,z)dydz = 0. (27) 

In Appendix 1X1 we prove the ergodicity of the process {z,y}, together with the 
fact that Co satisfies the Fredholm alternative. 

We will assume that the average of the velocity with respect to the invariant 
density p vanishes: 

(v(z)) p = 0. (28) 

From the identity J Td j Rd yC^p{y, z) dydz = and after an integration by parts 
using l|2t)|) . it follows that condition l|28|l is equivalent to 

(y) P = o. 

We assume that the following ansatz for the solution u e holds: 

u e = it + em + £ 2 u 2 + ... (29) 

with Ui = Ui(x, y, z,t), i — 1,2,.... We substitute l|29|l into l|24() and obtain the 
following sequence of equations: 



C uo = 0, 

CqUx = -Cxuo, (30) 

dug 

dt 



C u 2 = + 



From the first equation in l|3U|) we deduce that uo = uo(x,t), since the null 
space of Co consist s o f functions which are constants in y and z. Now the 
second equation in l|3l)|l becomes: 

Cou x = — l —y ■ X7 x ii . 



2 This pair z, y is the same as z, y solving 1231 up to a rescaling in time which is irrelevant 
to the ergodicity discussion here. 
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The centering condition that we have imposed on the vector field v(y, z) implies 
that (y) p = 0. Hence the above equation is well-posed. We solve it using 
separation of variables: 

ui = $(y,z) ■ VzUo 

with 

£ $(y,z) = -^y. (31) 
V T 

This is the cell problem which is posed on T d x R d . Now we proceed with the 
third equation in (|3*U|) . We apply the solvability condition to obtain: 

du 

= {L\ui) p 



at 



1 , , . d 2 u 



v r dxidxj 
= -^={y®<5>) p :D 2 x u . 

This is the Backward Kolmogorov equation which governs the dynamics on large 
scales. We write it in the form 

duo = K d 2 u 
dt 13 dxidx-j 



where the effective diffusivity is 



Jdj = A=(Vi*j) P . (33) 



IT 



The calculation of the effective diffusivity requires the solution of the cell prob- 
lem (|3ip. Notice that the cell problem is not ellipt ic - it is, however, hypo- 
elliptic. This follows from the calculations in Lemma [A. II 

Before studying various properties of the effective diffusivity, let us briefly 
present the b asic ingredients of the rigorous proof of the homogenization theorem 
presented in ^j- The basic idea is to apply Ito formula to the solution Q(y, z) 
of the cell problem to obtain, from (|23JI , 

1 f* 

x(t) = x{Q) H = / y(s)ds 



n 



(T 



x(0) - e ($(y(t), z(t)) - $(y(0), z(0))) + -= / V,$ (y(s), z(s)) dj3(s) 



' T 







It is straightforward to show that the bracketed terms on the right hand side of 
the above equation converge to 0, as e — > 0. By using the fact that the (y, z) 
process is fast, the martingale central limit theorem facilitates proof that the 
stochastic integral which appears in the above equation converges to a Brown- 
ian motion whose covar ian ce is given by the limit of the quadratic variation of 
the stochastic integral ^H], j^l Thm. 7.1.4]. The application of this theorem 
automatically provides us with the well-posedness of the limiting Backward Kol- 
mogorov equation. In the next section we establish this well-posedness directly, 
through the multiple-scales framework. 
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4 Properties of the Effective Diffusivity Tensor 

4.1 Well-posedness of the Limiting Equation 

In the previous section we showed that the dynamics of the inertial particles 
at large scales is governed by the backward Kolmogorov equation i|32|) . In this 
subsection we prove that this equation is well-posed, i.e. that the effective 
diffusivity is nonnegative. To see this, we first observe that 

2 2 

£*o(fp) = -p£of h — p\f + — V a p • Vyf, 

T T 

for every periodic function f(y, z) which is sufficiently smooth. We use this 
equation, together with an integration by parts and some algebra to obtain: 

2 f f 2 

f(C f)pdydz = -£-// \V y f\ 2 pdydz = ~?-{\V y f\ 2 ) p {M) 

lR d JT d ZT JR d JT d ZT 

This is the analogue of l|17(l for passive tracers. Now let cj> = a- $ where $ is the 
solution of the cell problem 1)3 ljl and a is a constant vector in R d . The scalar 
quantity <f> satisfies the equation 

£ o = — i=a-y. (35) 

From the formula for the effective diffusivity, together with equation (|34Jl . we 
obtain: 

a.K.a = —= [ J (a ■ y)(a ■ &)pdydz 

V T JM d Jf d 

4>£o4>pdydz 

'm d Ji a 

2t 



2 



V y (f>\ pdydz 

R d JT d 



(T 



2 



■^{\vA )p > °- 

4.2 Alternative Representation of the Effective Diffusivity 

The aim of this subsection is the derivation of an alternative representation for 
the effective diffusivity along the direction of the vector a in R d . To this end, 
we define the field \ thro ugh <p = y^ry.a + \ with (f> — <i>.a, $ being the solution 
of the cell problem (|31|l . Substituting this expression in 1)31)1 we obtain the 
following modified cell problem: 

C x = -a-v. (36) 

The effective diffusivity along the direction of the vector a expressed in terms 
of x i s: 

2 2 2 

aXa - °— \a\ 2 + ^-(\V yX \ 2 ) P + ^(^yX ■ a) p (37) 

Equ at ion s 1)36)1 and 1)37)1 have the same structure as the corresponding equations 
1)18)1 . I)2U)) for the first order dynamics. We will exploit this in the next section 
when we consider the r — > limit of 1)37) 1 . 
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4.3 Incompressible Flows 



We are unable to prove the analogue of what is known for passive tracers, namely 
that diffusion is always enhanced for incompressible flow fields. Numerical ev- 
idence, however, suggests that this enhancement is seen in a wide variety of 
situations and that, furthermore, the presence of inertia further enhances the 
diffusivity over the passive tracer enhancement. See section [S] 



4.4 Potential Flows 

From the representation 1)3 7|) we can prove that, for potential flows, the diffu- 
sivity is depleted for all t > 0, as is true for the case t = of passive tracers. 
As for passive tracers we use the fact that the expl icit form of the invariant 
measure is known for potential flows. From equation (|34|l . with / = 0, the facts 

that (\a.y\ 2 ) p = ^-|a| 2 and that a.\7 y p = —-^py.a and an integration by parts 
we obtain: 



u 



2 



= (V y x ■a) p =a-Ka - —\a 



2 



We use the above formula in equation l|37(l to deduce that, for potential flows, 
we have (a.\7 y x) P = — ^(|V a x| 2 )p which implies: 

a-K,a=^\a\ 2 -^-{\V yX \ 2 ) p . (38) 

Hence, transport is always depleted for potential flows. This formula should 
be compared with the formula (1221 arising in the case r = 0, passive tracers. 
We also remark that a more sophisticated analysis, based on the variational 
formulation of the effective diffusivity for passive tracers, yields that for potential 
flows, and at least in one dim ensio n, the diffusivity for r > is depleted even 
beyond its depletion for r = j3Ul Thm. 5.1]: 

9 

/C(r) < /C(r = 0) < °— . 

This is a sharper upper bound on /C(r) than the one that follows from (|38|) . 
On the other hand, equation (|38[) has the advantage that it provides us with an 
explicit expression for the differe nce b etween molecular diffusivity and IC(t) in 
terms of the solution of equation l|3tj[l . 



4.5 Conditions for the Centering Hypothesis (|2S|) to be 
satisfied 

In this sec tion we present some conditions which ensure that the centering con- 
dition It is easy to check that this condition is satisfied for potential flows. 
Indeed, the explicit form of the invariant density is known in this case and this 
enables us to perform the following computation: 

f v(z)p(y,z)dydz = —-r^rs I I W(z) e -^(^ 2 +^)) dydz 

fd J R d Zj {it (J*) 2 J t d J ffi d 

a 2 1 1 



2 Z (ira 2 ) 2 Jfd 7 K d 
0. 



V 



( e -Mh 2 +v(^ dydz 
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t 



Figure 3: Second moments of the particle position versus time for the velocity 
field l|4U|) . The graphs of t and t 2 are also included for comparison. 



In the general case, for which the invariant density is not explicitly known, we 
have to study the symmetry properties of the drift term in order to identify other 
classes of flows for which the centering condition is satisfied. Let us consider 
the case of a parity invariant flow, i.e. a flow satisfying the condition 

v (-z) = -v{z). (39) 

It follows from l|39|l that the solution of equation (|26|) -i.e. the invariant density- 
satisfies 

p(y,z) = p(-y,-z). 

Hence is satisfied. 

As a concrete exa mple which shows that condition (|39|l is necessary for the 
centering hypothesis to be satisfied, let us consider the two-dimensional 
velocity field 

v{z\, Z2) — (sin z\, — cos z±). (40) 

We have that 

vi(-z{) = -ui(zi) and v 2 (-zi) = v 2 (zi), 

and hence we expect diffusive long time behavior along the Z\ direction and bal- 
listic motion along the z 2 direction. In Figure 3 we present the second moments 
(xi(t) 2 }, (x2(t) 2 } of solutions to equation QJ with v(x) given by (f^Ufl 3 - 

3 The moments are obtained through Monte Carlo simulations. The details of the numerical 
simulations presented in this paper are discussed in section Iffl 
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4.6 Homogenization When the Centering Hypothesis (J5HJ) 
is not Satisfied 

In the previous section we imposed the centering condition l|28|) in order to en- 
sure that there is no mean drift and that the motion of the inertial particles at 
long scales is diffusive. In the case where this condition is not satisfied, then 
we expect that the effective behavior of the particles is described by a trans- 
port equation and that the diffusion appears only as a higher order correction. 
Indeed, an analysis similar to the one presented in the previous section, using 
the advective rescaling t = t/e and x = x/e, shows that in this case the parti- 
cle motion at large scales is goverened by the following backward Kolmogorov 
equation 

du du d 2 u 

dt 1 p dxi 13 dxidxj ' 

where (v) p = J Td J Rd v(z)p(y, z) dydz. Alternatively, the behavior of the parti- 
cles is diffusive at the reference frame moving with the mean flow. 



4.7 The Homogenization Problem for Equation (HJ) 

Let us now consider the equation Q in the presence of molecular diffusion. For 
steady flows this equation becomes 

x = Sv(x) ■ Vv(x) + - (v(x) - x + affj , (41) 

w here <5 is the ratio of fluid density to particle density. The techniques developed 
in |25| enable us to conclude that there exists a unique, smooth, invariant density 
for the process {x(t), ^/rx(t)}, which we denote by p s (y, z). It is straightforward 
to check that the assumption v(—z) = —v(z) ensures that the drift term b s {z) — 
Sv(z)-Vv(z) + ^v(z) is centered with respect to this invariant density, (b s ) p s = 0. 
Hence, the multiscale techniques developed in section [3] apply with —v(z) — > 

b s (z) and we can conclude that the rescaled process x e = ex(t/e 2 ), with x(t) 
being the solution of equation i|41|) . converges as e — ► to a Brownian motion. 
The covariance-effective diffusivity- of this Brownian motion is 

Vt 

where $ 5 solves the cell problem 

1 



C s 9 s = ~y, (42) 



with 



C s = ( V^b s ■ V v + -^y ■ V x ) + -C 



ou- 

T J T 



Now the effective diffusivity depends on 5 as well as on r and a. We study the 
effect of varying 5 in section by means of numerical experiments. The main 
interest here stems from the fact that, for a = 0, equation Q exhibits effective 
diffusive behaviour for certain cellular flows and choices of 8. How this diffusive 
mechanism interacts with molecular diffusion is a matter of some interest. 
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5 Small r Asymptotics 

In this section we show that the effective diffusivity tensor K. = K{t) reduces to 
the effective diffusivity tensor from first order dynamics, as r — > 0. To this end 
we define 

A =C OU , Ai = yV z + v(z)-V v 

and then, by polari zatio n, the effec tive diffusion tensor IC(t) is determined by 
l|37l) with x solving p solving i|26*|) and 

1 „ 1 „ 
Co = -A a + -=Ai. 

T V T 

We wish to expand x an d p in p owers of r and show that the leading order 
behaviour of JC{t) is given by (|20(l . Higher order terms in the small r expansion 
of the effective diffusivity JC(t) will be c omputed only for the case of a one 
dimensional potential flow in section \b. 41 Pertur bation c alculations similar to 
the ones presented in this section were reported in |431 144| for linear shear flows. 



5.1 Expansion for \ 

We set _ 

X = Xo + V^Xi + T X2 H 

in (EfuTl and then find 

Aoxo = (43) 
Aoxi = -AiXo (44) 
A0X2 = -a-v — AiXi- (45) 

The first of these equations implies that xo = Xo( z ) om y an d the second is 
soluble because 

{AiXa) p ou = (y ■ V z xo(z)) p ou = 0. 

Here p is the mean zero Gaussian invariant density of the OD process, satis- 
fying H'25f) . Solving for xi gives 

Xi = V ■ V zXo(z) + ipi{z) 
and the solvability condition for the xi equation yields 

YA z xo + v(z) ■ W z xo = -a.v(z). 

Thus xo is th e solution of the cell problem arising in the passive tracer case, 
equation (|18|l . Furthermore V y Xi — ^zXo( z ) an< i hence 

V yX = V^V zX o(z) + 0(T). (46) 

Notice that the fun ctio n ipi( z ) is undetermined at this point, but that it does 
not enter equation H46fl . 
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5.2 Expansion For p 

Now it remains to expand p(y, z) from 12(j|) in r. Notice that 

A* Q p ou (y) = 0, X7 yP ou (y) = ~2yp ou (y)/a 2 . 

Then, by with z —> y and Co — > Aq we find that 

Mf P OU ) = -P° U [A f ~ a 2 A,/ + 2y ■ V y f] 
= P OU Aof. 

Let p = p 0,7 7r. Then 



= L* (P OU *) 
= ±A* (p° U n) + ±Al(p OU tr) 

T y/T 

p ou p ou ( 2 

= Aqk H =■ y ■ V z 7r + v ■ X7 y n -(v ■ y)n 

t yjr \ er z 

Since p ou > everywhere we may divide through by it in the above expres- 
sion. If we then set _ 

TT = 7To + \ZriTl + T7T2 

we find that 

-Aqtv = 

2 

-Aottx = y ■ V^ttq + v ■ V a 7r -(v y)ir Q 

2 

-A ir 2 = y ■ V 2 7T! + v ■ WyTix - — (v ■ y)ir 1 . 

The first equation shows that 7To = tto(z) only and the second then gives 

2 

ni(y,z) = y ■ {V z tt 5-ttou} - pi. 

a 1 

As with the second term in the dual expansion, the function p\ is undetermined 
at this point. However, as the subsequent calculations will show, it is not needed 
in order to compute the first term in the small r expansion of the effective 
diffusivity. 

The negative of the right hand side of the third equation is 

2 2 
-y®y: {D 2 z ir a ? V z {tt v}} - v ■ {V z 7r ~ir v} 

2 , 2v r ^ 2 
+y®y : {v z n jn v\ <g> -~ + y ■ {V z p\ ?vpi\. 



Now 



a 2 



{y®y) p ou = —I and (y) p ou = 
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and so the solvability condition for W2 yields 

-yA 2 7r o + V 2 -{7r o t;} = 

and hence shows that ttq coincides with the passive tracer case, equation 
In summary, we have shown that 

p(y,z) = p ou (y)7T (z) + O(V^) (47) 

whe re tto(z) is the invariant density from the first order dynamics, satisfying 
CD- 

5.3 Limit of the Diffusivity Tensor 
Combining gSJl and gZJ in JSZjl gives 

2 2 p p 

a-Ka= ^-\a\ 2 + ^- / tt(z)\V zX o(z)\ 2 dz + a 2 ir(z)a.V zXq (z) + <D(s/¥) 
i i Jfd Jfd 

which, to leading order in r, is the expression for the first order dynamics. 

5.4 One dimensional Potential Flows 

The calculation of higher order terms in the small r expansion for the effective 
diffusivity is quite involved. Moreover, in the general case, these higher order 
terms do not seem to be of definite sign. However, it is possible to compute 
explicitly the next term in the small t expansion and to prove that it has a 
definite sign in the case of one dimensional potential flows and we now pursue 
this. Consider the equation 

tz = -V'{z) - z + a/3. (48) 

In this case we only need to solve perturbatively the c ell p roblem (|3 1 1) since 
the stationary Fokker-Planck equation corresponding to is exactly solvable 
yielding an invariant density independent of r: 

P (y,z) = ±^=e-Mh 2 + V(*)), (49 ) 

with Z = f Q exp(—-^V(zj) dz and y — -j^z. The effective diffusivity-which 
now is a scalar- is given by the formula: 

K=— / ycj)(y,z)p(y,z)dydz. (50) 

V r J-oo JO 

A lengthy calculation enables us to compute the first four terms in the small r 
expansion for the cell problem: 

<t>a = X(z), (51a) 

My,z) = y(x'(z) + i), (51b) 
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<h{v,*) = ly 2 x"(z) + Mz), (5ic) 

MV, *) = lv 3 x'"(z) + y (^"(z) - V'(z) X "(z) + ^(*)) . (51d) 

The cell pro blem for the first order dynamics, equation can be solved 

explicitly gf>| to give 



X(z) = / p(z) dz- z + c , 
Jo 

with 

p(z)=e^\ Z = t e^Uz. 



We use the above formula for x( z ) m equation (|15|l to obtain the following 
formula for the effective diffusivity 

nr = 0) = ^X. (52) 

The term tp2(z) satisfies an equation similar to the cell problem of the first order 
dynamics. The solution of this equation is 

M*) = fyv'iz)) 2 + P(z) dz + c x J' p(z) dz + ca, 

where 

c i = "4 l Z \{V'{z)fp{z)dz. 
v Jo 2 

The values of the constants Co, C2 are n ot n eeded for the computation of the 
effective diffusivity. We substitute now l|51l) and (|49|l . using the formulas for 
x(z) and ip2{z) and for the moments of the Ornstein-Uhlenbeck process. The 
final result is 

JC = £i-r f\{V'{z)fp{z)dz + 0{r 2 ). 
z ZZ Jo L 

The above formula shows that, for r sufficiently small, the diffusivity is depleted 
beyond the deplet ion exhibited by homogenization in the passive tracer case, 
which is given by 



6 Numerical Experiments 

In this section we study the dependence of the effective diffusivity for equations 
CJ, © on the non-dimensional parameters of the problem r, a and 5. For 
simplicity all the experiments we perform are for the Taylor-Green flow 

U(x) = ^ipTG jx), iPtg(x) = sin(xi) sin(x 2 ). 

4 In our derivation of the homogenized equation and the formula for the effective diffusivity 
we assumed that the velocity field is 1— periodic, rather than 27r-periodic. Of course the 
analysis, as well as the formulas that we derived, are trivially extended to encompass this 
change of period. 
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The closed streamlines of Lagrangian particle paths in this velocity field is a 
rather special situation and we describe numerical experiments for other stream 
functions, including open streamline topologies, in |34j . 

It is straightforward to check that the Taylor-Green flow satisfies condition 
l|39|) and hence the absence of ballistic motion at long scales is ensured. More- 
over, the symmetry properties of J5J imply that the two diagonal components of 
the effective diffusivity are equal, whereas the off-diagonal components vanish. 
In the figures presented below we use the nota tion K := JCn = K.22- 

Rather than solving the cell problem (|42[l. we compute the effective diffu- 
sivity using Monte Carlo simulations: we solve the equations of motion l|41|l 
numerically for different realizations of the noise and we compute the effective 
diffusivity through the formula 

K = lim ±-{{x{t) (x(t))) ® (x(t) - (x(t)))), 

t— >oo It 

where (•) denotes ensemble average. We solve the stochastic equations of motion 
using Milstein's method, appropriately modified for the second order SDE 19 
p. 386]: 

x n +2 = (2 - r)x n+ i - (1 - r)x n + rAtv(x n+ i) + crrAt7V(0, 1). 

where r = — and si = 1 — |, S2 = 1 + §.. This method has strong order 
of convergence 1.0. 5 We use N — 1024 uniformly distributed particles in 
27rT 2 with zero initial velocities and we integrate over a very long time interval 
(which is chosen to depend upon the parameters of the problem) with At m 
5.10- 4 min{l,r}. 

In some instances we compare the effective diffusivities for inertial particles 
with those for passive tracers. The latter are computed by solving the cell 
problem directly, by means of a spectral method similar to that described in 
|24| , together with extrapolation into parameter regimes where the dependence 
of the diffusivity is provably linear. 

6.1 The Effect of r on Diffusivity 

We compute the effective diffusivity as a function of the non-dimensional par- 
ticle relaxation time r for the Taylor-Green flow when 5 = 0.0. Our results 
are presented in Figure 0] for a = 0.1. For comparison, the effective diffusiv- 

2 

ity of the tracer particle (r = 0) and that of the free particle, namely are 
also plotted. The parameter t, apart from influencing the effective diffusivity, 
introduces an additional time scale into the problem |20j . In particular, for r 
large, we need to integrate the equations of motion over a longer time interval 
in order to compute accurate statistics. 

The main interest in this data is that it shows highly non-trivial dependence 
of the effective diffusivity on the parameter r as well as giving quantative infor- 
mation about how inertia enhances the diffusivity over that obtained in the case 
of passive tracers. We remark, in particular, that the effective diffusivity reaches 
its maximum for r = 0(1). It is in this regime that the free-flight time-scale for 
the inertial particle is of the same order as that induced by the velocity field. A 
sim ilar phenomenon, in a different context, has been observed by Vassilicos et. 
al. HSJ. 

5 All experiments reported in this paper have been independently verified by use of an 
alternative, linearly implicit, method. The agreement between the statistics computed using 
these two methods is excellent. 
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Figure 4: Effective diffusivity versus t for 5 = 0.0, a = O.f . 



6.2 The Effect of 5 on Diffusivity 

It is a well documented result E3 E3 EH1 EE] that the particle trajectories of 
(0J perform an effective Brownian motion even in the absence of noise, in certain 
parameter regimes. The linear stability analysis of (0J for the Taylor-Green flow 
indicates that in the parameter regime 5 £ 1) we expect a very complicated, 
chaotic behavior which might be interpreted as an effective Brownian motion 
at long times. When U(x) = V^^tgC^) then the noise free dynamics arising 
from when a = has two sets of equilibria, Xq = [nir, 0, mir, 0], n, m £ N 
and Xq = [§ + nn, 0,5+ mir, 0], n, m £ N. The first set of equilibria have two 
dimensional unstable manifold for 6 > r _1 and one dimensional for S < t ; 
the second set have a two dimensional unstable manifold for S < 1 and is stable 
for 5 > 1. Numerically we observe diffusive behaviour, in the absence of noise, 
if and only if S £ (t , 1). In this regime, all equilibria have two dimensional 
unstable manifolds. The basic mechanism for diffusion is chaotic mixing caused 
by separation near the separatrices of equilibria; it is unsurprising therefore that 
the linear stability of equilibria play a strong role on determining the interval 
of 5 in which diffusion occurs. 

In Figure^ we plot the second moment of the particle position as a function 
of time for values of 5 below and above the threshold S = 0.5, and in the absence 
of noise. As expected, the particle motion is ballistic for S — 0.45 and diffusive 
for S = 0.55. 

In Figure Eh we plot the effective diffusivity as a function of 5, again in 
the absence of noise. Since we have chosen r = 2.0 we expect an effective 
Brownian motion for <5 6 (0.5, 1). Notice that the effective diffusivity increases 
as S — > 0.5 + , and appears to diverge in the limit. This is to be expected, since 



22 



10' 



11 



10 




Figure 5: Effective diffusivity as a function of S for 
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8 

Figure 6: Effective diffusivity as a function of 6 for a = 0.1 and r = 2.0. 




a. 5 = 0.0, t= 1.0 b. 5 = 0.7, r = 2.0 

Figure 7: Effective Diffusivity versus a for 5 = 0.0 and (5 = 0.7 . 
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S = 0.5 separates ballistic motion (for 6 < 0.5) and diffusive motion (for S > 0.5) 
motion. 

In Figure |H1 we plot the effective diffusivity as a function of S when the 
particles are subject to additional molecular diffusion. In this case the effective 
diffusivity is a decreasing function of 5. 

6.3 The Effect of a on Diffusivity 

It is well known, see e.g. (22], that for the case of passive tracers the effective 
diffusivity depends on the molecular diffusivity a in a highly non-linear, very 
complicated way. In particular, the limit as a tends to is singular and the en- 
hancement in the diffusivity-for divergence free flows- depends crucially on the 
topology of the streamlines. It is therefore interesting to study the dependence 
of the effective diffusivity on a for the inertial particles problem. 

In Figure we present the effective diffusivity as a function of the molecular 
diffusion for two sets of parameters: (a) for S = 0.0 (for which the streamlines 
are closed), r = 1.0 and (b) for S — 0.7, r = 2.0 (a regime in which there exists 
a well defined effective diffusivity even in the respectively. The diffusivity of the 

free particle is also plotted for comparison. Moreover, in Figure[7|[a), we also 
plot the effective diffusivity for passive tracers. 

In both figures we see the clear enhancement of diffusivity over the bare 
molecular value. Furthermore, in Figure[7{a) we also observe that the enhance- 
ment in the diffusivity is significantly greater for r > (inertial particles) than 
it is for t = (passive tracers), especially when a is small. 

A further interesting observation concerns the dependence of the diffusivity 
on a. For <5 = the dependence is highly non-trivial, exhibiting both a local 
maximum and a local minimum. For S = 0.7 the presence of inertia leads to 
an effective diffusivity which increases as a becomes smaller. This should be 
contrasted with the established fact that, for passive tracers in periodic cellular 
flows, t he effective diffusivity decreases linearly in er, for a sufficiently small 
jnniE3- It would be interesting to understand how the enhancement scales with 
d — ► and to compare it with known theoretical results for the passive tracer 
case; in particular, it would be in tere sting to extend the theory of maximally 
and minimally enhanced diffusion |24) to the inertial situation studied here. 

7 Conclusions 

The problem of periodic homogenization for inertial particles is considered in 
this paper. It is shown that, at long times and large scales, the inertial particles 
perform Brownian motion and a formula for the effective diffusivity is derived. 
Furthermore, the dependence of the effective diffusivity on the non-dimensional 
particle mass r a(Stokes number) and the ratio of the fluid to particle density 
S is studied, by means of analysis and numerics. 

It is shown that, as r — » 0, the effective diffusivity converges to the one 
obtained from the homogenization of passive tracers. Moreover, it is shown 
through numerical experiments that for a variety of interesting divergence free 
flows, the diffusivity in the presence of inertia is enhanced much beyond the well 
documented enhancement of the diffusivity for passive tracers. Furthermore, the 
dependence of the effective diffusivity on r and 5 is studied numerically in some 
detail. 

The calculation of t he eff ectiv e diffusion tensor requires the numerical so- 
lution of equations l|3T|l and It is coneivable that this task might be as 
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computationally demanding as direct Monte Carlo simulations, because the do- 
main of the PDE is unbounded in the momentum variable, and because the 
PDE is not elliptic (only hypo-elliptic). This is to be contrasted to the case of 
passive tracers in periodic flows; there th e ca lculation of the effective diffusivity 
requires the solution of the elliptic PDE on a periodic domain. Equations 
of this type can be routinely and effeciently solved using, for example, a spectral 
method. From this point of view our results might not provide any computa- 
tional advantage over Monte Carlo simulations. However, the results reported 
in this paper provide a mathematical framework for rigorous analysis of the 
dependence of the effective diffusion coefficient on the physical parameters of 
the problem. We have already undertaken such an analysis to study the limit of 
small Stokes number and we plan to undertake further studies in future work. 

The numerical results reported in section for a simple two- dimensional 
steady flow, exhibit a wide range of interesting physical phenomena. As ex- 
amples we mention the dependence of the effective diffusivity on the Stokes 
number and the fluid/particle density ratio, for a given streamline topology, 
question which are, we believe, of great interest to the applied community. The 
purpose of this work is to develop a framework within which such questions 
can be addressed. We plan to investigate some of these issues in future work. 
The dependence of the effective diffusivity for a wider class of velocity fields is 
undertaken in |34|. 

Summarizing, we note the following specific areas where future work would 
be of interest: 

• the extension to time dependent velocity fields v(x,t), either periodic in 
time or ra ndom in tim e - for example with an Ornstein-Uhlcnbcck struc- 
ture as in |3~3"1 14TH |4T| : 

• the extension to random velocity fields in space; 

• rigorous analysis of the parametric dependence of the effective diffusivity 
on <7, r and <5, taking into account the free streamline topologies; 

• further numerical studies for velocity fields other than the simple Taylor 
Green flows studied here - in particular to study problems where the 
Lagrangian particle paths have open streamline topologies, such as the 
Childress-Soward family; this is initiated in |34| . 



In this appendix we prove the existence of a unique invariant measure for the 
process {x,x} solving (JTJ and, moreover, that the generator C of the process 
satisfies the Fredholm alternative. This justifies the formal multi-scale calcula- 
tions presented in section [3] To simplify the notation we set a = r = 1. The 
equations of motion become 
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Appendix 



x = v(x) — x + 0. 



(53) 



The generator of the Markov process {x, y\, with y = x is 



C = yX7 x + ( v (x) -y)-V y + -A. 



■y 



We have the following theorem 
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Theorem A.l. Assume that v(x) £ C°°(T d ). Then there exists a unique, 
smooth invariant density p(x,y) for the process {x, y}: 

£*p(x,y) = 0. 

Let further h{x, y) be a smooth function such that J„ d J Td h{x, y)p(x, y) dxdy = 
0. Then the Poisson equation 

-Cf = h (54) 

has a unique mean zero solution in L 2 (T d x M. d ,e~ s dxdy) for every 5 6 
(0,2a- 2 ). 

The proof of the exi sten ce of a unique invariant measure for our process is 
based on the results of [25] and is broken into three lemmas. First, we need to 
prove the existence of a smooth transition probability density for our Markov 
process. This is accomplished by means of Hormander's theorem Thm 
V38.16]. Then we need to prove the compactness of phase space, for which we 
need to find an appropriate Lyapunov functions. Finally, we need to prove that 
the transition probability density is everywhere positive. To show this we need 
to use a controllability argument. Th e proof of the existence and uniqueness of 
solutions of the Poisson equation (|54|l is based on Fredholm's theorem. 

Lemma A.l. The Markov process generated by C has a smooth transition prob- 
ability density. 

Proof. This follows by an application of Hormander's theorem. The basic idea 
behind this theorem is that, even though noise does not act directly to the 
position variable, there is nevertheless sufficient interaction between momentum 
and position so that noise, and consequently smoothness, is transmitted to all 
degrees of freedom. We write the generator in Hormander's "sum of squares" 
form: 

1 - 

i=l 

where Xi = i = 1, . . . d and X — y - V x + (v(x) — y) . Let now [A, B] denote 

the commutator between the vector fields A, B and let Lie{F} denote the Lie 
algebra generated by the family of vector fields F. Define 

A = Lie{Xi, . . .X d } 

and 

A k = Lie{LY , U], U e A-i}, k = 1, 2, . . . 

Set finally 

H = Lie{^o,^i---}- 

According to Hormander's theorem, a sufficient condition for the Markov process 
{x, y} to possess a smooth invariant density is for 7i to span the tangent space 
T x ^ y M 7 where M = T d x R d . We readily check now that 

d d 

[X ,X t } = -— + —, i = l,...d, 
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and consequently 

Span (Lie {A Q , A%}) = T x , y M. 

Thus, Hormander's hypothesis is satisfied and the Markov process generated by 
C has a smooth density. □ 
Now we prove that the existence of a Lyapunov function. 

Lemma A. 2. There exists a constant (3 > such that the function V(x,y) = 
1 + ^11 2/11 2 satisfies 

C(V(x,y))<~V(x,y)+f3. 

Proof. We have that V(x, y) maps the state space onto [1, oo) and that limjjyjj^oo V (x, y) 
oo. Moreover we have 

C(V(x,y)) = v-y-\\y\\ 2 + ^ 

1 „ „ 2 ( d 1 l|2 



< - 2 \\vr + [- 2 + - 2 \M 

< -V(x,y)+[3, (55) 

with (3 = d/2 + (l/2)sup xeTd ||w(.t)|| 2 + 1. □ 
The last ingredient which is needed for the proof of the ergodicity of the 
process generated by C is the fact the transition probability P t is everywhere 
positive. 

Lemma A.3. For all z := (x, y) £ T d x R d , t > and open O C T d x R d , the 
transition kernel for (|53Jl satisfies Pt(z,0) > 0. 

For the proof of t his le mma we refer to \2bl Lemma 3.4]. 

Proof of Theorem ^. 1\ The existence of a unique invariant mea sure follows 
from Lemmas I A. 1 1 1 A .21 and I A . 3l upon using Corollary 2.8 from In order to 
prove the existence and uniqueness of solutions of the Poisson equation l|54|) we 
ne ed t o prove that the generator C has compact resolvent. This is accomplished 
in [H Thm 3.1, Thm 3.2]. □ 

Remark A.l. The above lemmas enable us to conc lud e that the system, con- 
verges exponentially fast to its invariant distribution ]2S$ . 

Remark A. 2. We a bit of extra work we can also prove sharp estimates for the 
invariant distribution and the solution of the Poisson equation —Cf = h. We 
refer to fJ 61 Thm 3.1, Thm 3.2] for details. 
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